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Abstract 

We numerically investigate the formation of D-dimensional black holes in high-energy particle 
collision with the impact parameter and evaluate the total cross section of the black hole produc- 
tion. We find that the formation of an apparent horizon occurs when the distance between the 
colliding particles is less than 1.5 times the effective gravitational radius of each particles. Our nu- 
merical result indicates that although both the one-dimensional hoop and the (D — 3)-dimensional 
volume corresponding to the typical scale of the system give a fairly good condition for the hori- 
zon formation in the higher-dimensional gravity, the {D — 3)-dimensional volume provide a better 
condition to judge the existence of the horizon. 
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I. INTRODUCTION 



The brane world scenario is paid much attention in the context of the unified theory 
of elementary particles. This scenario regards our space as the 3-brane with large extra 
dimensions on which gauge particles and interactions are confined. In this scenario, the 
Planck energy can be at the O(TeV) scale [|I|. One of the consequences of lowering the 
Planck scale is that the properties of a small black hole whose radius is smaller than the 
size of the extra dimensions are substantially altered; the black holes is well-described as 
a D-dimensional black hole centered on the brane, but extending out into the {D — 4)- 
dimensional extra space and the radius of such a black hole is ~ 10'^^ times larger than that 
of the usual black hole with the same mass. Hence it becomes much easier to produce black 
holes using a future planned accelerator such as the CERN Large Hadron Collider and this 
possibility has been discussed by many authors 0. 

The possibility of producing black holes by the collision of particles can be estimated 
by using the hoop conjecture 0, which states that an apparent horizon forms when and 
only when the mass M of the system gets compacted into a region whose circumference C 
satisfies 

(J 

^ WM) - ^' 

where rh{M) is the Schwarzschild horizon radius for the mass M. If we assume that the 
inequality (|I]) give the condition for black hole formation in the higher- dimensional gravity, 
we can evaluate the impact parameter h of the colliding particles which leads to black 
hole production. The circumference that surrounds two particles at the instant of collision 
is C ~ 2h. By setting 2b/2TTrh{2fi) ~ 1 where 2/i is the center of mass energy of the 
system, the maximal impact parameter femax that leads to the black hole formation becomes 
^max ~ ''"/i(2yu). By introducing a numerical factor F{D) close to unity, the total cross section 
for black hole production is written as the following form: 

CTb.h.production = F {D)nrl{2fi) . (2) 

Obtaining F{D) is necessary to improve experimental predictions of the black hole produc- 
tion in collider physics and observations of ultra-high energy cosmic rays. 

In our previous paper 0] , we investigated the black hole formation in high-energy head-on 
collisions of particles in the Z)-dimensional gravity and found that Hd becomes a parameter 



2 



to judge the existence of the apparent horizon for all D. Although this implies that the hoop 
conjecture holds for this system, we do not know the condition for the horizon formation 
in the higher-dimensional spacetime. Recently, Ida and Nakao |^ investigated momentarily 
static initial data in the five- dimensional spacetime and found that a spindle distribution 
of matter can lead to the black hole formation even if the hoop C is arbitrarily larger than 
27rr/i(M). This provides a counter example for "only when" part of the hoop conjecture. 
They showed that the isoperimetric inequality (characteristic area) < 47rr^(M) is satisfied 
on the horizon and conjectured that for a characteristic {D — 3)-dimensional volume Vd_3 
of the system, the inequality 

would become a condition for the horizon formation in the D-dimensional spacetime with 
the gravitational constant Go- This "(Z^ — 3)-volume conjecture" (volume conjecture) is 
one candidate that provides a better condition for the horizon formation in the higher- 
dimensional spacetime than the hoop conjecture and it is worth investigating. 

In this paper, we consider the grazing collisions of particles in the Z)- dimensional Einstein 
gravity and investigate the formation of apparent horizons. Eardley and Giddings devel- 
oped a method of finding apparent horizons for this system. The problem was reduced to 
a boundary-value problem for a Poisson's equation in a (D — 2)-dimensional flat space and 
they solved it analytically for D = 4 case. We solve this problem numerically for D > 4 and 
obtain the maximal impact parameter fomax for black hole formation and the factor F{D) in 
Eq. (Ip. Then using these solutions, we discuss the condition of horizon formation for this 
system from the view point of the hoop conjecture and the volume conjecture. 

This paper is organized as follows. In Section II, we briefly review the method of finding 
apparent horizons in the system of the grazing collision of particles. In Section III, we 
present our numerical results and discuss the two conjectures for this system. Section IV is 
devoted to summary and discussion. 

II. APPARENT HORIZONS IN THE GRAZING COLLISION OF HIGH- 
ENERGY PARTICLES 

To simplify the situation, we consider black holes with horizon radius smaller than the 
size of extra dimensions. This enables us to ignore the effect of the brane tension and the 
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geometry of the extra dimensions. The metric with a high-energy point particle is obtained 
by infinitely boosting a Schwarzschild black hole with fixed total energy fi. The result- 
ing system becomes a massless point particle accompanied by a plane-fronted gravitational 
shock wave which is the Lorentz-contracted longitudinal gravitational field of the particle. 
Combining two shock waves, we can set up the collision of high-energy two particles moving 
in ±z direction. This system was originally developed by D'Eath and Payne |^ and recently 
analyzed in @, P]. The metric of this system outside the future light cone of the colliding 
shocks is given as 



HtP = + ^e(n)V.V,<l>(a. - 0.+), (4) 

where u = t — v = t ^ G is the Heviside step function and x = (x*) is the point in 
fiat [D — 2)-space (xi, ■ ■ ■ , xd-2) that is transverse to the direction of particle motion. The 
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function (^{x) depends only on r = \x\ = \fx^ and takes the form 

$(a;) = -8G4yulogr, for D = 4 (5) 

where fi_D-3 is the volume of a unit [D — 3)-sphere and x± denote the points of two particles 
in [D — 2)-dimensional space and we take 

x^ = (±6/2,0,- ■■ ,0). (7) 

The apparent horizon S is defined as a closed spacelike {D — 2)-surface on which the 
outer null geodesic congruence has zero convergence. Eardley and Giddings developed a 
method of finding the apparent horizon in the union of the two shock waves, u = > v and 
V = > u. Their method reduces the problem to finding the {D — 3)-dimensional closed 
surface C and two functions ^±(a;) on the {D — 2)-dimensional plane satisfying 

V2 _ = interior to C, (8) 
^± = on C, (9) 
• =4 on C, (10) 
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where $± = $(a; — x±). If we can find a solution of Eqs. (|),® and (|T0|), a surface with 
zero expansion is given as the union of two {D — 2)-surfaces v = — \&+(a;) in m = > f and 
u = — \&_(a3) in f = > M. C is the intersection of the two {D — 2)-surfaces at m = f = 0. 
For D = 4, Eardley and Giddings solved the above problem analytically They used the 
conformal invariance of two-dimensional Poisson's equation. For D > 4, we cannot use their 
method because (D — 2)-dimensional Poisson's equation does not have conformal invariance. 
We must solve the above problem numerically. 

In the (D — 2)-dimensional space, C and have a symmetry about xi-axis that connects 
two source points. Hence solving the apparent horizon is reduced to a two-dimensional 
problem. We use spherical coordinate (r, 6) to represent the point x in {D — 2) -dimensional 
space, where r is the distance |a;| from the origin and 6 is the angle between x and Xi- 
axis. We parameterize the curve C hj r = g{9). The function g{6) and have symmetry 
g{6) = g{7i — 6) and ^+(r, 6) = \E'_(r, tt — 6). We calculate C and numerically as follows. 
First, we give a trial curve C and solve the equation (|D under the boundary condition @. 
Then, we calculate the value 6 = ■ V^- — 4 on C, and determine a modified boundary C 
using the value of 6. C and \E'+ converge to the solution of Eqs. (j^) (j^) (0) by iterating these 
two steps. In order to proceed the first step, we introduce a regular function h = — ^+ 
that obeys D-dimensional Laplace equation 

h,rr + — — -hr- + + ^ cot 6hg = 0. (11) 

The boundary condition is 

= 0, (12) 



n\ = cE> I — 

"'\r=g{e) QQ 



_ dh 
e=o 06 



that comes from (^) and the symmetry about xi-axis. To simplify the boundary condition, 
we introduce a new coordinate (f, 9) by 

r = g{e)f, 6 = 6. (13) 

In this coordinate, the equation (|1T]) becomes 

g'i J '^"^ r2'^,ee ~ ^ ".ef 

1 / q" q'^ ~q'\ D-A 

+ - [ D - ?> - — + 2^ - iD ~ A) coi6- \ h f + — ^ cot 6h ^ = Q (14) 
^ V 9 r 9 J ' 
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and the boundary C is given by r = 1. We solve (|T^) using the finite differential method 
with (50 X 100) grids. There is a coordinate singularity and we cannot write down the 



where p is the distance from xi-axis. We relate the value at the grid points around f = 
using Eq. (p!5|). 

Once h converges, we proceed the second step to determine the modified boundary C. 
We calculate 6{9) = V\l/i ■ V\l/2 — 4, and determine the modified boundary C as 



If e is sufficiently small, C converges. In our calculation, we took the value of e as typically 
10^'^ ~ 10^^ times the horizon radius of the system. Smaller e is required to converge 
the calculation for the larger impact parameter b. The maximal impact parameter 6max is 
obtained when dg{6,b)/db becomes — oo. For D = 10 and 11, the numerical instability 
occurs at the neighborhood of 6 = femax- We obtained fomax by extrapolating g{Tc/2,b) by 
fitting the numerical data. 

To evaluate numerical errors, we compare the numerical solution with the analytic solu- 
tion given by Eardley and Giddings for D = A case. The numerical error is less than 0.02%. 
For D > 5, we evaluate the error by comparing the solution with the result of the calculation 
with (100 X 200) grids. The error is about 0.04% for D = Q, 0.07% for D = 8, 0.25% for 
D = 10, and 0.4% for D = 11. 

III. NUMERICAL RESULTS AND CONDITION FOR HORIZON FORMATION 

Figure 1 shows the shape of C on (xi, X2)-plane for Z) = 4, 7. For 6 = 0, the apparent 
horizon consists of two {D — 2)-dimensional fiat disks whose radius tq is 



The radius g{6) of C takes the maximum value rmax at 6' = 0, vr and takes the minimum value 
Tniin at ^ = 7r/2. As the impact parameter b increases, the radius of the {D — 3)-surface C 




r = g{9) + ex 6{9). 



(16) 




(17) 
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FIG. 1: The shape of the apparent horizon C on (xi, X2)-plane in the colhsion plane u = v = ^ 
for D = 4,..., 7. Incoming particles are located on the horizontal line X2 = 0. Values of b/rQ are 
0(»), 0.4(o), 0.7(4), 1.0(0, shown only in D = 6,7), and 6max/'^o(*)- As the distance b between 
two particles increases, the radius of C decreases. 

becomes smaller. The shape of C at 6 = bmax deviates from a sphere and the ratio rmax/'^mm 
becomes large as D increases and more prolate shaped horizon can exist for larger D. 

Figure 2 shows the relation between b and rmin for each D. The value of b^g^^/ro ranges 
between 0.8 and 1.3 and becomes large as D increases. Because rh{2fi) ^ tq, 6max/'''/i(2/i) is 
also 0(1) and increases with D. This behavior of &max/'"h(2/i) is different from the case of 
the head-on collision [Q, where the ratio (distance of two particles)/r/j(2yu) at the horizon 
formation decreases with D. Our numerical results for 6max/2ro and F{D) are summarized 
in Table I. 

Now we investigate the condition for horizon formation in the two particle system. Be- 
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FIG. 2: The relation between the impact parameter b and the minimum radius rmin of C for 
D = 4, 11. The value of bjaax/i"o grows as D increases. 



TABLE I: The value of 6max/2ro and F{D) for D = 4, 11. 



D 


4 


5 


6 


7 


8 


9 


10 


11 


Vax/2ro 


0.402 


0.480 


0.526 


0.559 


0.583 


0.603 


0.619 


0.632 


F{D) 


0.647 


1.084 


1.341 


1.515 


1.642 


1.741 


1.819 


1.883 



cause 6max/''"/i(2/i) ~ 0(1), the hoop conjecture can be apphcable to estimate the maximal 
impact parameter in high-energy colhsions in the higher-dimensional spacetime. But the 
hoop conjecture cannot explain the increase of &max/^/i(2/i) with the increase of D. Figure 
3 is the plot of b^g^^^D). We found that the relation b^nax — l-5r/i(/i) well describes the 
behavior of bjaax{D) obtained by numerical calculation. As r/i(/i) is an effective gravitational 
radius of each incoming particle, the horizon formation occurs when there is an overlap of 
more than one half of r/i(/i) between the regions with radius r/i(/i) around each particles. 
Because r/i(/i) is proportional to the ratio rh{^)/rh{2^) = 2^-*^/^^"^^ increases with 

D. This leads to the increase of &max/''"/i(2yu) with D. 

The behavior &max/'''/i(2/u) oc 2-1/(^-3) can be derived from the more fundamental prin- 
ciple, the [D — 3)-dimensional volume conjecture. In the present system, the condition 
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FIG. 3: The value of 

^max (crosses) as a function of the spacetime dimension D. The dotted Une 
rft(2/i) which comes from the hoop conjecture. Although the ratio ^max 

/rfi{2fi) takes the value 

around unity, the hoop conjecture does not explain the increase of 6max/''h(2/x) with D. The solid 
line bma.x/rh{2fj,) = 1.5r/,(^)/r/j(2^) ~ 2"i/(^~3) gives a good fit of 6max- 



can be written as 



n 



Vn- 



l/(D-3) 



This condition reduces to the hoop conjecture (0) for D = 4 case. As the characteristic 
length of the system is b in xi-direction and rft(/i) in other directions, the characteristic 
{D — 3)-voIume becomes Vds ~ brj^~'^{ij) and Hd becomes ~ [brj^~^{jj)/rj^~^{2fi^~^^^''^ 



By setting TYd ~ 1, the maximal impact parameter is given by 

D-4 



(19) 



r„(2^) 

Although this is rough order estimation, the calculation indicates that the volume conjecture 
provides a better condition than the hoop conjecture. 

To test the validity of the two conjectures further, we consider how much energy is 
trapped by the black hole. This is related to the definition of the mass in the system. In 
our previous paper [Q, we discussed the existence of the gravitational wave energy that 
does not contribute to the horizon formation in the system with motion and found that the 
hoop conjecture with Hawking's quasilocal mass provides a good condition for the horizon 
formation in four-dimensional gravity. In the present system, we calculate the quantity 



Ma.h. = 



D 



n 



D-2 



(20) 
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FIG. 4: The relation between the horizon mass Ma__iiJ2ii and the impact parameter 6 for D = 
4, 11. MA.H./Syu becomes small as D increases. The black hole in the higher dimensional 
spacetime can trap small amount of the energy. 

where A0-2 is {D — 2) -dimensional area of the apparent horizon. By the area theorem, 
-^A.H. provides the lower bound of the final mass of the black hole and coincides Hawking's 
quasilocal mass on the horizon. Therefore Ma.h. can become an indicator of the energy 
trapped by the horizon. Figure 4 shows the behavior of Ma.h./S/x as a function of b for each 
D. We find that the value of Mp^n /2fi at 6 = femax decreases as D increases. In the higher- 
dimensional spacetime, the amount of "junk" energy increases because the gravitational field 
distributes in the space of the extra-dimensions and only a small portion of the total energy 
of the system can contribute to the horizon formation. This junk energy will be radiated 
away rapidly after the formation of the black hole. This strongly suggests that we should 
use the quasilocal mass to check the hoop and the volume conjectures. 

To test the validity of the two conjectures with the quasilocal mass, we must calcu- 
late M(S), C(S) and Vd-3(S) for all surfaces S, and then take the minimum value of 
C(S)/27rr^(M(S)) and ¥0-3(8) /nD-3rh~^{M{S)). In the present system, calculation of the 
quasilocal mass is difficult because we treat the black hole on two connected null hyperplanes, 
while ordinary quasilocal mass is defined in spacelike hypersurface. Hence we calculate the 



value of f/^-H- = C(SA.H.)/27rr;,(MA.H.) and = [Vd^S a.r.) /^Dsr^-' (Ma.h.)] 



on the apparent horizon Sa.h.- In calculating \ we must specify the hoop for the 
surface Sa.h.- The hoop is defined as the longest closed geodesic on the surface. In our 
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FIG. 5: The value of H^-^- function of 6/r/j(2^) for D = A, 11. The circles show the values 
at 6 = 6max- ^^'^'(^max) bccomes large as D increases. 
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FIG. 6: The value of H^-^- function of 6/r/j(2/i) for D = 4, 11. The circles show the values 
at 6 = 6max- '^^•^■(6max) becomes small as D increases. 

system, the hoop is the length of the curve C on (a;i, a;2)-plane shown in Figure 1. Figure 5 
shows the numerical results of H^'^ {b) for each D. H^'^ {b) monotonically increases with 
b for all D and it can be a parameter to judge the existence of the horizon. The value 
at horizon formation -ff^'^ (&max) ranges between 1.4 and 1.6 and increases with D. The 
deviation of the value Hd from unity is due to the definition of the mass and we expect 
that this value approaches close to unity if we use the appropriate mass of the system and 
the hoop conjecture holds. The increase of -ff^'^ (&max) with D reflects the fact that more 
prolate shaped horizon can form for larger D and 6max/'"h(2/i) increases with D. 

As for {D — 3)- volume, general definition has not existed yet. Ida and Nakao^ selected 
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some {D — 3)-surfaces by considering symmetries of the horizon and took the maximum 
value of their {D — 3)-volume. Using the similar method, we adopt the {D — 3)-volume of 
the {D — 3)-surface C. Figure 6 shows the numerical results of Ti^^'ih). 'H^-^ {b) is almost 
constant and slightly increases with b: 



(21) 



Thus the isperimetric inequality 'H^^'ip) < 1 holds on the horizon in this system. The 
factor in Eg . (^T]) contradicts the constancy of Ti;^'^' (^max) with respect to D that we have 
expected from (plQl). It may be absorbed to the definition of the volume factor, but we do 
not have the rigorous definition of the volume and we cannot explain the meaning of this 
factor at this stage. By ignoring this factor, the behavior of Ti^^ ib) for each D is similar 
to H^'^ {b) for D = 4. This indicates that Tin with the quasilocal mass provides the value 
that is close to unity at the horizon formation and the volume conjecture holds. 



IV. SUMMARY AND DISCUSSION 



We have investigated the black hole formation in the grazing collision of high-energy 
particles numerically. The black hole is produced when the impact parameter b satisfies the 
condition b < 1.5r/i(/i). This condition can be derived by using the volume conjecture and 
suggests that the volume conjecture provides a better condition for horizon formation than 
the hoop conjecture in the higher- dimensional gravity. 

We evaluated the value of and Tio on the apparent horizon and both inequalities 
Hd ^ 1 and T^d ^ 1 give fairly good conditions for the horizon formation. This result 
is in contrast to the horizon formation by the spindle matter distribution 0; In this case, 
the gravitational field in the transverse direction of spindle matter is {D — l)-dimensional 
and the black hole forms even if the hoop C is far greater than 27rr/j(M) ioi D > 5 and 
Hd ^ 1 does not give the condition for the horizon formation. But in the system with two 
point particles, the gravitational field in the middle of two particles is weak if the distance 
of two particles are larger than rh{M) and becomes sufficiently strong to form the horizon 
if the distance of two particles equals to rh{M) for any D. This leads to the condition of 
the black hole formation Ho ~ 1. We therefore expect that the hoop conjecture in the 
higher-dimensional spacetime holds for the system that consists of two point particles. 
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The behavior of T[^'^ {b^s.x) indicates that 'Hd with the quasilocal mass provides the 
better condition for the horizon formation compared to the hoop conjecture. But what we 
have confirmed is the necessary condition for the horizon formation; if an apparent horizon 
exists, there is a surface that satisfies Ti/j < 1. There remains the possibihty that even if the 
apparent horizon does not exist, there is a surface that satisfies the TY/j < 1. Hence further 
investigation is required to confirm whether {D — 3)-volume conjecture provides a sufficient 
condition for the horizon formation. This is our remaining problem. 
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